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ABSTRACT 


In  mechanizing  a computer  to  solve  the  lead  equations  for  an  airborne 
flexible-gun  fire-oontrol  system,  accurate  values  are  needed  for  the 
input  quantities  of  range  and  sight-line  angles-off  and  their  derivatives. 

It  is  the  purpose  of  this  report  to  investigate  the  effectiveness  of  various 
eleotrioal  filters  and  the  theories  underlying  them  in  simultaneously 
smoothing  and  differentiating  various  input  signals  provided  by  the  radar. 
Emphasis  has  been  placed  upon  evaluating  RC-networks. 


CONCLUSIONS 


A partial  investigation  of  the  effectiveness  of  several  eleotrioal 
filters  and  the  theories  underlying  them  has  been  made . In  partioular , 
it  has  been  found  that,  for  smoothing  and  differentiating  the  angular 
coordinates  of  the  sight-line  as  provided  by  the  traoking  radar,  a filter 
based  upon  the  Zadeh-Ragazzini  theory  is  no  more  effective  than  a single 
stage  RC-filter  in  the  case  of  linearly  varying  time  signals  and  noise 
characterized  by  an  exponentially  deoaying  autocorrelation  funotion.  This 
oonolusion  is  still  valid  for  non-linear  time  signals  if  the  single  stage 
RC-filter  is  replaoed  by  an  RC-network  characterized  by  a transfer  function 

s(l  + a.s  + aas*  + ...  + a ,sn"1) 

H (s) i — T , 

1 + a,s  + a^s2  + ...  + a s“  +bs  + ...  + b,sJ 
1 n-x  n j 

with  j > n.  Beoause  of  this  fact,  the  main  effort  in  this  report  deals 
with  transfer  functions  of  this  type. 

The  investigation  has  been  partial  in  that  the  signal  input  functions 
were  of  two  types  only  and  the  noise  autocorrelation  funotion  a deoaying 
exponential.  However,  the  two  signals  were  representative  of  the  two  most 
probable  target  paths:  pursuit  courses  and  straight  line  interception 

attaoks.  The  noise  autocorrelation  function,  selected  on  the  basis  of 
available  radar  traoking  data,  is  fairly  representative,  muoh  more  so  than 
the  assumption  of  white  noise  implicit  in  the  Blaokman-Bode-Shannon  theory. 
For  this  reason  the  latter  theory  has  not  been  investigated  further. 

In  general,  it  appears  that  RC-networks  with  transfer  functions  H L, (s) 

nj 

oan  be  used  to  give  results  within  required  limits  for  major  portions  of 
nearly  all  realistic  tactical  courses.  The  design  of  a range  rate  filter, 
considered  briefly  in  this  report,  must  await  definite  knowledge  concerning 
radar  range  noise  characteristics.  The  problem  here,  however,  is  eased  by 
the  fact  that  higher  order  derivatives  of  range  rate  are  small  and  as  a 
result  signal  distortion  oan  be  kept  within  bounds  more  easily  than  in  the 
oase  of  an  angular  rate  filter. 


ii 


CONFIDENTIAL 


CONFIDENTIAL 


SECURITY  INFORMATION 


MAVQRD  REPORT  1424 


TABU!  OF  CONTENTS 

FOREWORD i 

ABSTRACT  . ii 

CONCLUSIONS ii 

I.  INTRODUCTION 1 

A.  STATEMENT  OF  THE  PROBIEM 1 

B.  SIGNALS  USED 1 

C.  NOISE  ASSUMPTIONS  3 

D.  CRITERIA  FOR  FILTER  EFFECTIVENESS 4 

II.  SOME  PARTICULAR  THEORIES  OF  FILTERING 6 

A.  GENERAL  BACKGROUND 6 

B.  ZADBH-RAGAZZINI  THEORY 7 

C.  BLACKMAN,  BODE  AND  SHANNON  THEORY 9 

D.  THE  RC  DIFFERENTIATING  FILTER 9 

III*  EVALUATION  OF  RC  ANGULAR-RATE  FILTERS 13 

IV.  EVALUATION  OF  RC  RANGE- RATE  FILTERS 25 

NOTATIONS 29 

REFERENCES 31 

FIGURES 

1.  GUN-TARGET  SPACE  AND  RELATIVE  PATHS 2 


CONFIDENTIAL 


SECURITY  INFORMATION 


iii 


TABUS 


1 . COMPARISON  07  Z-R  AND  RC  A-FII/TBRS  12 

2.  A-FILTBR  - OUTPUT  ERRORS  - H,,(s)  STRAIGHT  LINES  CASK 

(p  - 90°) 11 . 16 

3.  A-FILTKR  - OUTPUT  ERRORS  - H,  '(s)  STRAIGHT  LINES  CASE 

(p  - 90® ) 38 18 

4.  A-FILTER  - OUTPUT  ERRORS  - H,-(s)  STRAICHJT  LINES  CASE 

(p  « 90°) 20 

5.  A- FILTER  - OUTPUT  ERRORS  - Max.  A0  Min.  Q STRAIGHT  LINES 

CASE  (p  * 135  ) . . . ' 21 

6.  A-FILTER  - OUTPUT  ERRORS  -/ Na  - 100  mil8a/seo8  lfox.  A0 

o 

STRAIGHT  LINES  CASE  (p  m 135°) 22 

7.  A-FIl/TBR  - OUTPUT  ERRORS  - Na  = 25  milsa/seoa  Max.  A0 

0 

STRAIGHT  LINES  CASE  (p  = 135°) 23 

8.  A-FILTER  - OUTPUT  ERRORS  - H_,(s)  STRAIGHT  LINES  CASE 

(P«135°* . . 24 

9.  r-FILTBR  - NOISE  OUTPUT 26 

10.  r-FILTER  - OUTPUT  ERRORS  - H_,(s)  STRAIGHT  LINES  CASE 

(P  = 135°) ?3 27 


CONFIDENTIAL 


confidential 


NAVORD  REPORT  1424 


I . INTRODUCTION 


A.  STATEMENT  OF  THK  PROBIBM. 

In  mechanizing  a oomputer  to  solve  the  lead  equations  for  an  air- 
borne flexible-gun  fire-oontrol  system,  aoourate  values  are  needed  for 
the  input  quantities  of  range  and  sight  line  angle-off  and  their 
derivatives.  The  range  r is  furnished  direotly  by  the  radar  while 
the  sight  line  angle-off  is  approximated  by  the  azimuth  and  elevation 
angles,  A and  E,  respectively,  of  the  axis  of  the  radar  antenna.  The 
quantities  r,  A and  E thus  furnished  are  not  smoothly  varying,  being 
contaminated  with  noise.  The  noise  is  generated  as  the  sum  of  several 
different  effects  which  may  be  noted.  These  are 

a.  Wandering  of  the  radar  beam  over  the  target  surface, 

b.  Fading  of  the  radar  echo  as  a result  of  ohanging  target  aspeot 
and  atmospheric  variations , 

c.  Antenna  and  reoeiver  noise, 

d.  Servo  follow-up  noise. 

The  derivative  quantities  r,  A and  B,  obtained  by  differentiating  the 
rough  r,  A and  E values,  will  depart  sharply  from  the  desired  derivative 
values  unless  steps  are  taken  to  smooth  the  original  input  data. 

It  is  the  purpose  of  this  report  to  investigate  the  effectiveness  of 
various  electrical  filters  and  the  theories  underlying  them  in  simultan- 
eously smoothing  and  differentiating  the  azimuth  angle  A of  the  axis  of 
the  radar  antenna  as  it  tracks  a target  moving  on  a prescribed  path 
relative  to  the  ownship.  The  quantity,  A(t),  free  of  noise,  will  be 
spoken  of  as  the  signal.  The  noise,  superimposed  upon  the  signal,  will 
be  denoted  by  N(t).  *" 

B.  SIGNALS  USED. 

Two  different  types  of  time  functions  were  used  for  A(t),  one 
associated  with  the  target  on  a straight  line  path  relative  to  the  own- 
ship  and  the  other  associated  with  the  target  on  a pure  pursuit  course 
relative  to  straight  line  motion  of  the  ownship.  For  the  first  of  these, 
the  spaoe  path  of  the  target  was  assumed  to  form  an  angle  p with  the 
path  of  the  ownship,  as  shown  in  Figure  1. 
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FIQDRS  1.  GUN-TARGBT  SPACE  AMD  RKLATITB  STRAIGHT  LIMB  PATHS 

Denoting  the  ownship  and  target  speeds  by  and  VT>  the  relative 

speed  by  Y,  the  cross-over  range  by  rQ  and  rQ/V  by  k,  we  note  that 

the  funotional  form  of  the  signal  for  azimuth  angle , measuring  time 
from  oross-over,  is 

(1)  A(t)  « G - tan"1  | 

The  time  interval  of  interest  for  the  rear  turret  of  a bomber  would 
be  that  preceding  cross-over.  The  interval  considered  is  -8  < t < 0. 
Data  used  were  YQ  = 175  yd/sec.,  V T = 225  yd/seo.,  = 400  yardsT  with 

C and  k being  determined  by  combining  these  with  p . 
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Initial  evaluations  of  filters  vrere  made  assuming  p ■ 90°.  Later 
evaluations  vere  based  upon  a compromise  assumption, of  p m 135°. 

The  seoond  signal  used  assumed  a parabolic  fit  to  A(t)  for  A(t) 
associated  with  the  target  on  a pure  pursuit  course.  In  particular, 
with  units  in  radians  and  seconds , 

(2)  A(t)  = .0125 ta  - .It  - .2  (0  < t < 8) 

The  maximum  angular  rate  here  amounts  to  200  mils/seo.  at  t ,■  4. 

C.  NOISE  ASSUMPTIONS . 

The  noise  function  N(t),  associated  with  A(t),  is  not  oapable  of 
direct  analytic  formulation  and  henoe  must  be  dealt  vith  statistically 
using  the  basic  concepts  of  autocorrelation  and  spectral  density.  (Tor 
a good  introduction  to  these  techniques,  the  reader  should  oonsult 
Chap.  VI  of  reference  (1).)  From  numerous  souroes  (see  references  (1), 

(2) , (3))  a very  typical  form  of  the  autocorrelation  function  R(t)  for 
noise  in  angle  A,  obtained  experimentally  in  extensive  radaur  tracking 
tests, is 

(3)  R(t)  « a*  e"altl  . 


The  quantity  <j?  = R(o)  gives  the  mean-square  value  of  the  noise 
amplitude  N(t), 

T 

(4)  a*  ■ lim  ^ J Na(t)dt. 

-T 


This  follows  direotly  from  the  "time-average’'  definition  of  R(t), 
namely 


(5)  R(t) 


. lim 
T 


L 2T  / 

-T 


N(t)  N(t  + t)  dT 


The  spectral  density  funotion  G(f)  for  the  input  noise  N(t)  is, 
by  the  Wiener-Khinohine  theorem,  the  Fourier  oosine  transform  of  R(t). 
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The  funotion  0(f)  multiplied  by  df  gives  a direct  measure  of  the  amount 
of  noise  power  in  the  frequenoy  range  f 
tation  (3)  for  R(t),  ve  find  that 


(6)  0(f)  - 4 / R(t)  oos  2nft  dt 

o 

where 

(0  as  Snf. 

The  oonstants  a and  a3  appearing  in  (3)  depend  upon  the  radar 
type,  the  rate  of  change  of  target  aspeot  and  other  factors.  Values 

of  a in  the  referenoes  oited  show  variation  from  5 to  40  (seo  1). 

In  this  report  a * 5 and  a = 10  are  aooepted  as  representative 
values.  For  a3  the  value  7.1  (mils)3  was  taken  in  agreement  with 
the  value  given  in  reference  (2/. 

D.  CRITERIA  FOR  FILTER  EFFECTIVENESS . 

A general  outline  will  now  be  given  of  the  criteria  used  in  deter- 
mining the  effectiveness  of  a given  filter  in  smoothing  and  differen- 
tiating the  input  A(t)  + N(t).  Considering  the  noise  N(t)  with 
spectral  density  0(f),  we  desire  to  compute  the  mean-square  value  of 

the  output  noise,  whioh  we  designate  by  the  symbol  N3  • (The  sub- 
script o indicates  "output",  the  dot  refers  to  the  faot  that  the 
input  noise  is  differentiated,  and  the  bar  yields  the  mean  of  the 
quantity  N3  .)  If  we  denote  by  GQ(f ) the  spectral  density  of  the 
output  noise  N0(t),  then  by  the  faot  that  fi3  = Rq(0),  where  RQ(t) 
is  the  autocorrelation  function  for  NQ(t) , and  the  Viener-Khinohine 
relation, 

RQ(t)  = / GQ(f)  oos  2nft  df  , 


we  find  the  all-important  relation 


(7)  Nj  - f Go(f ) df  . 


to  f + df.  With  the  represen- 


4a  a3 
a3  + m3 
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To  find  GQ(f)  we  need  the  theorem  relating  GQ(f  ),  G(f)  and  the  filter 
transfer  function  H(s).  [See  reference  (1),  p.  288.]  This  is 

(8)  G0(f)  = G(f)|H(s)|8  for  s = j©  . 

The  transfer  function  H(s)  is  in  turn  defined  to  be  the  ratio  of  the 
Laplaoe  transforms  of  output  and  input  functions  to  the  filter.  An 
equivalent  definition  describes  the  transfer  function  H(s)  as  the 
Laplace  transform  of  the  filter  weighting  function,  W(t).  Thus, 


H(s)  = r W(t)  e"st  dt  . 
o 


The  funotion  H(t),  sometimes  known  as  the  smoothing  function  or  memory 
function,  is  the  response  of  the  filter  to  the  unit  impulse  function 
6 (t ) as  input.  The  latter,  known  also  as  the  Dirad  delta  funotion,  is 
defined  by  the  relations 


m 

6(t)  * 0 for  t 0 , 6(0)  = oo  , / 6(t)dt  » 1 . 


Combining  (6)  and  (8)  we  have 

(9)  N*  - / G(f)|H(j©)|a  df  . 


Consider  now  the  signal  input  to  the  filter,  A(t).  The  desired 
filter  output  is  A(t),  the  actual  output  is  AQ(t).  We  shall  denote 

the  signal  distortion,  Aq  - A by  the  symbol  A0 . The  quantity  A6 

will  be  a funotion  of  t and  of  filter  parameters  Fx , Fa , etc. 

(For  a simple  capacitance-resistance  filter,  e.g.,  Fx  ® RC#)  For 
reasons  indicated  later,  all  the  F's  are  assumed  equal.  Thus  we 
write: 

(10)  Signal  Distortion  = A0(F,t). 

It  now  remains  to  combine  (9)  and  (10)  to  yield  a criterion  for  filter 
effectiveness.  To  this  end,  the  following  two  definitions  are  made. 
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(11)  A9m  » Uaz.  Value  of  A0(F,t)  with  respect  to  t 


— — n A 

(12)  A0a  - i / A0a(F,t)  dt  = Mean  Square  Value  of  A0(F,t) 


(T  is  the  time  on  the  path,  here  taken  to  be  8 seoonds.) 
Relations  (9),  (11)  and  (12)  are  now  oombined  as  follows: 

(13)  Q - / A0j  + Nj 


(14)  Q = / A6a  + Na 


It  should  be  noted  that  A0„,  A02  and  Na  are  each  functions  of  P alone. 

M ' O 

On  the  basis  of  the  two  oriteria  (13)  and  (14),  filters  are  compared  by 

noting  how  small  their  corresponding  Q's  or  Q' s become  for  a given  type 

of  signal  input.  The  criterion  using  O’ is  somewhat  more  realistio  than 
that  using  Q. 


H . SOUS  PARTICULAR  THEORIES  OP  FILTERING 


A.  GENERAL  BACKGROUND. 

Basing  their  work  upon  the  fundamental  work  of  Wiener  (see  refer- 
ence (4)),  numerous  writers  have  put  forth  theories  of  filtering  where- 
in the  filter  is  characterized  by  a weighting  function  which  is  the  best 
possible  according  to  some  definite  mathematical  criterion  for  selecting 
the  "best"  funotion  out  of  a whole  olass  of  funotions.  In  particular 
there  should  be  mentioned  the  fundamental  papers  by  Phillips  and  Weiss 
(reference  (5)),  Cunningham  and  Hynd  (reference  (6)),  and  Zadeh-Ragazzini 
(reference  (7)).  Sinoe  the  last  of  these  is  more  general,  including  the 
other  papers  as  speoial  oases,  we  shall  include  a very  brief  outline  of 
it. 
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6.  ZADBH-RAQAZZINI  THEORY. 

In  the  Z-R  theory,  the  signal  is  assumed  to  be  representable  as  a 
polynomial  P(t)  of  degree  not  higher  than  a speoified  number  n.  The 
noise  funotion  N(t)  is  assumed  to  be  a stationary  function  of  time 
described  by  the  autocorrelation  funotion  R(t).  [For  a preoise  defi- 
nition of  "stationarity" , see  reference  (l),  p.  270.]  Denoting  the 
weighting  funotion  of  the  filter  by  l(t),  where  l(t)  is  defined  for 
0 <.  t <,  T*  and  is  assumed  zero  elsewhere,  we  may  write  for  the  filter 
output,  ®0(t),  corresponding  to  the  input  E^Ct)  «*  P(t)  + N(t): 

T* 

(15)  B0(t)  = / W(T)[P(t  - t)  + N(t  - t)]  dx  . 

o 

If  the  moments  of  1(t) , 


(16)  « / xr  l(x)  dx  , r * 0,1,2,. . . ,n  , 


are  introduced  and  P(t  - x)  is  expanded  as  a polynomial  in  t,  we  may 
rewrite  (15)  in  the  form: 


(17)  Ko(t) . ^op(t)  - ^p(t)  + f(t)  + ...  + (-i>n  ^ p(n)(t) 

T* 

+ / l(x)  N(t  - x)  dx  . 
o 

Consider  now  the  error  of  the  filter,  namely, 

e m Output  - Desired  Output; 
or  in  symbols, 

(18)  e(t)  - Bo(t)  - P(t) 


The  Z-R  theory  now  seleots  a "best  possible"  filter  as  one  whioh  is 
characterized  by  a W(t)  determined  to  satisfy  the  following  two 
oonditions : 
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(•)  li*  r / e(t)  dt  * 0 , 

L-m  U 0 

L 

(b)  a0*  *=  lim  i / e2(t)  dt  is  a minimum. 

L—m  " o 


Noting  (17)  and  (18),  we  may  write  condition  (a)  in  the  equivalent  form: 


(19)  P(t)  - fi0P(t)  - nxP(t)  +|‘P(t)  + ...  + (-if  ^ P^t)  . 


(Equation  (19)  assumes  that  the  noise  N(t)  has  zero  mean.) 

Equation  (19)  is  equivalent  to  imposing  n + 1 constraints  on  the 
function  W(t): 

= 0,  i «■  0,2,3,...  ,n;  (i  / 1) 

(20) 

Hi  = -1  • 

It  is  of  interest  to  rewrite  (18)  in  the  light  of  (20). 

T* 

(21)  e(t)  * / l(x)  N(t  - t)  dr 


Since  the  filter #is  linear  and  differentiating,  we  reoognize  the  right 
side  of  (21)  as  NQ(t),  and  henoe  ».(t)  « NQ(t).  Condition  (b)  is  thus 

equivalent  to  minimizing  RQ(0),  the  mean  square  value  of  the  output 

noise.  Thus,  « RQ(0). 

For  the  aotual  minimization  of  RQ(0),  the  condition  (b)  can  be 
rewritten  as 

Ji  J# 

(22)  Rq(0)  • / / l(t)  1(t)  R(t  - x)  dt  dt  . 

o o 
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The  Minimization  procedure  involves  an  application  of  the  Caloulus  of 
Variations.  The  result  will  now  be  given  for  the  autocorrelation 
assumption  made  in  (3).  It  is  found  that 

(23)  W(t)  > Aq  + Ajt  + A2t2  + ...  + Antn  + Cx6(t)  + Dx6(t  - T*), 

where  the  A' s , Cx , Dx  are  constants  which  depend  on  T* . The  acutal 

minimum  value  of  R (0)  = Na  for  the  case  n = 1 turns  out  to  be 

o'  7 o 

(24)  F « 24-  ■*  °* 

° T*(aaT*a+  6aT*  + 12) 


C.  BLACKMAN,  BODE  AND  SHANNON  THEORY. 

In  addition  to  the  work  cited  in  paragraph  A of  this  section,  there 
should  be  mentioned  the  paper  by  Blackman,  Bode  and  Shannon  (reference 
(8))  which  bases  the  optimum  filter  on  a weighting  function  W(t)  which 
minimizes  the  mean  square  prediction  error  of  the  filter  under  the 
assumption  of  "whits"  noise,  i.e.  noise  with  a speotral  density  which 
is  oonstant  for  all  frequencies  (flat  spectrum).  The  reason  for  this 
assumption,  according  to  B-B-S,  is  that  actual  noise  spectra  are  "sub- 
ject to  variations  due  to  faotors  which  it  is  not  desirable  in  praotioe 
to  attempt  to  control" . The  B-B-S  theory  oan  be  considered  to  be  a 
particular  oase  of  the  Z-R  theory. 


D.  THE  RC  DIFFERENTIATING)'  FILTER. 

It  is  possible  to  oompare  directly,  using  the  criterion  developed 
in  equation  (4),  the  effectiveness  of  filtering  achieved  by  a Z-R 
filter  and  that  achieved  by  a network  of  constant  resistances  and 
oapaoitanoes , a so-called  RC-f liter.  In  particular  it  oan  be  shown 
that  the  simple  RC-filter  shown  in  Figure  2 is,  under  fairly  reason- 
able assumptions  as  to  signal  and  noise,  entirely  equivalent  to  an 
n ■ 1,  Z-R  filter. 
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FIGURB  2 

To  see  this,  let  us  compare  the  minimum  values  of  Q for  the  two 
filters  (see  eq.  (13)  for  the  definition  of  Q).  The  expression  for 

N2  for  the  Z-R  filter  is  given  by  (24).  The  corresponding  expression 

for  the  RC-filter  is  found,  using  (9),  to  be 


(25)  *0  - ■ <7  * R0)- 

[ The  noise  is  assumed  to  be  desoribed  by  the  autocorrelation  function 
given  by  (3).) 

Now  for  the  Z-R  filter, 


»\a  24  a o8  a c3 

N — 11 ■■■■'■'  = --  ■ — -■  » ■ - ■—  ■■  ■ - 

0 T#  (12  + 6aT*  + a2T*a ) (^)[l  + a (^)  + &8^‘ j 

whioh,  for  small  T*  and  moderate  values  of  a,  is  approximately 

w m [-•(§]] 

mt 

Thus,  we  note  that  (26)  is  the  same  funotion  of  £ that  (25)  is  of  F. 

Let  us  now  oompare  the  signal  distortions  produced  by  the  two  filters, 
assuming  the  signal  to  be  non-linear  in  time.  We  find  for  the  RC-filter: 


10 


CONTIDBNTIAL 


0 


or,  upon  expanding  A(t  - t)  in  a Taylor  series  and  integrating, 

(27)  A9  = -F  A(t)  + Fa  A(t)  - F3  AU)(t)  + ... 

Since  F is  generally  small  and  A(t)  and  higher  derivatives  of  not 
much  importance,  ve  see  that 

(28)  A0  = -F  X(t)  . 

From  (17)  ve  note  that  for  n = 1,  Z-R, 

A0  = ^ A(t)  - ^ A(t)  + ...  . 

A direct  calculation  shows  that  jj^  *s  -T* , independent  of  a.  Hdnoe , for 
assumptions  already  made, 

(29)  A0  i - (§*]  X(t)  , 

T# 

which,  again  ve  note,  is  the  same  funotion  of  g that  (28)  is  of  F. 

Ve  are  thus  able  to  oonolude  that  the  value  of  F vhioh  minimizes 
Q(F)  is  the  same  as  the  g vhioh  minimizes  Q(T#). 

Computations  vere  made  using  the  signals  described  in  section  I.B. 
The  results  appear  in  Table  1 vith  representing  the  minimum  of  Q 
with  respeot  to  F or  T* . 
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III.  EVALUATION  OF  RC  ANGULAR-RATS  FILTERS 


In  the  prooess  of  investigating  the  synthes izati on  of  oirouits 
designed  to  approzinate  the  weighting  functions  for  Z-R  differentiating 
filters,  Washington  University  Researoh  Foundation  (see  referenoe  (9)) 
found  that  nearly  the  same  results  oould  be  achieved  with  proper 
combinations  of  resistances  and  oapaoitanoes  without  the  oomplioations 
encountered  in  inoluding  the  required  delay  lines. 

The  transfer  functions  of  the  differentiating  RC  networks  oan, 
in  general,  be  expressed  by 

s4l  axs  + a^s  + ...  + a.s  ) 

(30)  H .(s)  « - i r-  (n  - i + 1,  J > n) 

1 + b-^s  + b2s  + bgS  + ...  + b.s 

u 

The  emphasis,  in  evaluating  filters,  has  been  plaoed  upon  the 
formulation  of  the  transfer  funotion  and  its  effect  upon  the  toted 
errors  Q and  Q as  defined  by  (13)  and  (14),  respectively.  The  Laplaoe 
transform  of  the  normal  response,  Eo(t),  to  an  input  signal,  B^lt),  is 
given  by 


(31)  L{B0(t)}  « Hn  J(s)  L{B1(t)} 
or 

(32)  B0(t)  W(T)  B^t  - T)  dT 

where 

(33)  W(t)  - L"l(Hnti(s)}. 

If  H ,(s)  is  expanded,  "s"  treated  as  a differentiating  operator 

uj 

and  the  inverse  Laplaoe  transform  found  when  the  input  signal  is  A(t) 
and  the  output  AQ(t),  it  is  noted  that 

(34)  AQ(t)  « A(t)  - \k(t)  + kgA^t) 
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In  this  expression,  the  k are  determined,  by  the  network  constants 

r 

ar  and  b^.  One  notes  that  setting  a^  = br  for  r = l,2,3,...,i 
reduces  this  to 

(35)  AQ(t)  = i(t)  ~ b^+1jA^  + ^(1+2)^^  “ •••• 


where  the  coefficients  of  the  higher  order  terms  are  still  determined 
by  the  network  constants. 

Since  the  desired  output  is  A(t),  it  is  seen  that  the  distortion  is 
given  by 


(36) 


A0  . = - b 
nj 


(i+i) 


A(i+2)(t) 


+ C 


(i+2) 


A(i+3)(t) 


This  reduces  to  the  result  given  in  (27)  when  ax=  0,  bx  = F and 

a = b — 0 for  r = 2,3,...  . 
r r 


If  it  happens  that  the  input  signal  can  be  expressed  as  a polynomial 
of  degree  not  greater  than  n,  then  (36)  shows  that  the  signal  distortion 
will  vanish  identically.  Thus , it  is  seen  that  a network  whose  transfer 
function  is  given  by 


s (1  + b jS  + b 2s  2 + ...  + bn  xsn~  X) 


(57)  H^s). 


1 + bjS  + bjS* 


+ b.s^ 


(j  > n) 


will  provide  the  undistorted  derivative  of  a polynomial  signal  of  degree 
n or  less . 

,i  It  is  very  unlikely  that  input  signals  could  be  represented  entirely 
by  a single  polynomial  of  low  degree;  however,  it  is  noted  that,  if  the 
higher  order  derivatives  do  not  become  large,  the  distortion  can  be  made 
small  by  proper  adjustment  of  the  network  constants.  For  the  major 
portion  of  realistic  attaok  courses,  these  higher  order  terms  are  small; 
however,  for  interception  courses  (straight  lines)  the  higher  order  terms 
become  large  near  cross-over.  This  type  of  course  provides  the  supreme 
test  in  evaluating  a filter. 

After  the  distortion  is  determined,  it  must  be  oombined  with  noise, 
as  found  by  using  (9),  to  find  Q and  Q.  Adjusting  the  network  con- 
stants to  reduce  distortion  is  found  to  inorease  the  noise  output  and 
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The  tabular  results  indioate  that,  for  these  networks  and  this 
particular  signal,  the  total  error  for  most  points  would  be  far  too 
large;  however,  a comparison  of  the  results  for  the  straight  line  and 
pursuit  oourss  in  Table  1 indicates  that  the  total  error  can  be  reduced 
considerably  for  other  courses.  IJ  was  deoided  that  a compromise 
straight  line  oourse , with  p « 135° , would  be  more  realistio  in 
evaluating  networks  and  that  network  transfer  funotions  with  extra 
terms  in  the  denominator  should  be  used  to  reduoe  the  noise,  sinoe 
noise  was  the  major  contributor  to  the  total  error.  It  was  also 
decided  that  the  major  portion  of  future  investigations  of  this  typp 
of  network  should  be  done  experimentally,  sinoe  it  is  believed  to  be 
the  more  economical  approach. 

To  indicate  the  regions  to  be  investigated,  some  results  are  pre- 
sented wherein  the  previously  described  method  for  approximating  the 
maximum  distortion  was  used. 


TABUS  5 

A-FILTER  - OUTPUT  ERRORS  - Max.  AG  liin.  Q 
STRAIGHT  LINES  CASE  (P  = 135°) 


Transfer 

Function 

Pole  for 
Uin.  Q 

A9II 

mils/sec 

mils2/seca 

Q 

mils/seo 

H22 (s  ) 

-3.5 

10.5 

232.7 

18.5 

H23 (s ) 

-6.3 

9.7 

182.2 

16.6 

H24(s) 

-8.2 

11.4 

193.7 

18.0 

H33(s) 

-2.4 

8.7 

256.7 

18.2 

H34(s) 

-3.8 

8.75 

170.4 

15.7 

^35  (s  ) 

-5.1 

9.05 
- — 

216.3' 

17.3 

The  results  of  Table  5 indicate  that  Q still  does  not  approach 
the  3 mil/seo  total  error,  which,  at  present,  is  considered  as  an 
upper  limit.  It  may  be  found  that  Q can  be  brought  into  this 
neighborhood. 
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Sinoe  noise  is  always  present,  it  seems  more  logical  to  reduoe  its 
R1IS  value  to  near  the  upper  allowable  limit  on  the  total  error  and  then 
examine  the  distortions  for  a number  of  realistic  taotical  courses. 
Perhaps  the  A0's  for  a large  number  of  these  oourses  oan  be  oombined 
statistically  or  the  expeoted  distortion  oan  be  determined  from  a 
statistical  representation  of  all  realistic  oourses.  Tables  6 and  7 
present  the  maximum  distortion  and  total  error,  Q , to  be  expeoted  when 
the  equal  poles  are  fixed  so  as  to  make  noise  equal  to  100  milsa/seo* 
and  25  mils3/seea,  respectively.  A0U  was  computed,  as  before,  from 
the  expansion.  u 


TABUB  6 

A-FILTER  - OUTPUT  ERRORS  - N * = 100  mils2/seo* 

l 

Max.  A0  - STRAIGHT  LINES  CASE  (p  = 135°) 


Transfer 

Function 

Poles 

• 9 

No 

A0^ 

Q 

Haa(s) 

-2.14 

100 

28 

29.7 

Ha3(s) 

-4.75 

100 

17 

19.7 

Hm(s) 

-6.11 

100 

20.6 

22.9 

H33(s) 

-1.4 

100 

43,7 

44.8 

H3,(s) 

-3 

100 

17.8 

20.4 

H3S(s) 

-3.7 

100 

23,7 

. 25.7 

22 


CONFIDENTIAL 


CONFIDENTIAL 


HAVORD  REPORT  1424 


TABUS  7 

A-FILTER  - OUTPUT  ERRORS  - N*  = 25  mill? /sec? 
liax.  A0  - STRAIGHT  LINKS  CASE  (p  = 135°) 


Transf er 
Funotion 

Poles 

K 

*eu 

Q 

Haa(s> 

-1 

25 

128 

128.1 

Ha3(s) 

-2.64 

25 

55 

55.2 

Ha4(s) 

-3.5 

25 

62.7 

62.9 

H33(s) 

-0.66 

25 

417.4 

417.5 

H34(s) 

-1.7 

25 

97.7 

97.8 

HjbCs) 

-2.15 

25 

120.7 

120.8 

From  the  foregoing,  it  would  appear  that  the  networks  described  by 
transfer  functions  H23(s)  and  f^4(s)  would  give  better  results;  how- 
ever, the  ewidenoe  is  not  conclusive  since  this  analysis  is  based  upon 
necessary  approximations . 

As  information  of  possible  interest,  we  include  the  numerical  resulj 
for  two  R^a(s)  networks,  evaluated  for  the  straight  line  case,  p = 135 

with  unequal  poles  and  no  approximations  made  in  determining  signal 
distortion  from  (36). 
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TABUS  8 

A-FILTBR  - OUTPUT  ERRORS 
STRAIGHT  LINES  CASS  (R 


„ , \ s (81  + 58.5s) 

H*3<s>  * (s +"3 J (s' + 4 . 5 ) (s  +6) 

t 

46  mils/sec 

-8 

-.348 

-7.5 

-.427 

-7.0 

-.530 

-6.5 

-.663 

-6.0 

-.839 

-5.5 

-1.070 

i -5.0 

-1.376 

-1.780 

-2.310 

-3.5 

-2.984 

-3.0 

-3.784 

-2.5 

-4.585 

-2.0 

-5.010 

-1.5 

-4.229 

-1.0 

-.962 

-0.5 

5.601 

0 

13.658 

Na  mil*/seo8 
0 

80.8 

Q 

9.7 

Q 

16.6 

24 


A 


46  mils/sao 


-0.984 
-1.201 
-1.477 
-1.832 
-2.295 
-2.901 
-3.699 
-4.758 
-6.122 
-7.936 
-10.160 
-12.659 
-14.820 
-15 .118 
-10.787 
1.113 
20.001 


23.1 

9.6 

20.6 
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IV.  EVALUATION  OF  RC  RANGE- RATE  FILTERS 

In  analyzing  networks  designed  to  provide  smoothed  range  rate  data, 
the  emphasis  has  been  placed  on  network  transfer  functions  of  types 
Hgj(s).  Since  the  third  and  higher  order  time-derivatives  of  range 

are  small  for  reasonable  tactical  courses  of  interest  in  designing  a 
rear  turret  defensive  system,  it  was  thought  unnecessary,  and  probably 
not  profitable,  to  resort  to  more  complicated  networks.  For  a more 
extensive  investigation  along  these  lines  see  reference  (10). 

As  in  the  angular  rate  case,  higher  order  time  derivatives  of  range 
were  found  to  be  larger  for  interception  (straight  line)  oourses  than 
for  pursuit  type  oourses.  For  this  reason  it  was  deoided  to  plaoe  the 
emphasis  upon  the  compromise  straight  lines  course,  with  P » 135°.  If 
we  measure  time  from  cross-over,  we  have 


(38)  r = y/ 160,000  + 25,565  t2 


(-8  < t < 0) 


Since  almost  no  information  concerning  the  nature  of  noise  corrupting 
radar  range  data  was  available , emphasis  was  placed  upon  distortion  in 
the  output  signal.  It  was  desirable  that  the  total  error  in  range  rate 
be  less  than  5 yd/sec.  For  purposes  of  analysis,  it  was  deoided  that 
networks,  having  half  this  total  error  as  their  maximum  distortions, 
should  be  designed  and  tested  for  attenuation  of  an  assumed  noise,  These 
approximate  maximum  distortions  were  found  by  the  same  methods  used  in 
the  angular  rate  oase.  Transfer  functions  of  networks,  having  2.5  yd/seo 
maximum  distortion  with  the  prescribed  signal  input,  are  given  by 


(39) 


H22(s)  = 

s (8 .75  + 6s  ) 

F 

+ 2.5)(s  + 3.5) 

H23(s)  = 

s(l46 ,545  + 88.07s) 

F 

+ 3.5)(s  + 5 .3 )(s  + 7.9) 

H24(s)  = 

s (3003. 7392  + 1667.136s) 

F 

+ 5.4)(s  + 6 .6  )(s  + 8 .6 )(s  + 

9T87 


Temporarily,  an  assumption  was  made  for  noise  characteristics  to 
facilitate  a comparison  of  these  three  filters.  This  assumption  is 
based  on  the  findings  reported  in  reference  (12).  From  this  we  have 
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100  G. (0) 

(40)  G (f ) — i-  where  a>  « 2nf 

1 102  + a8 


The  above  reference  indicates  that  most  of  the  noise  in  ranging  is 
probably  due  to  "glint".  G^(0)  *=  kL3  for  glint  noise  where  k and  L 

are  dependent  upon  the  aspect  and  type  of  target.  A compromise  is  made 
which  seems  to  be  in  keeping  with  present  day  tactics.  The  attacking 
plane  is  assumed  to  be  near  head-on  in  aspect  and  both  planes  are  assumed 
to  be  doing  some  maneuvering.  We  have  assumed  G^ (0)  = 7 x 10_aL3  ft3/cps 

where  L is  the  length  of  the  attaoking  plane.  In  referring  to  measure- 
ments of  modern  fighter  planes , it  was  found  that  L may  be  as  great  as 
50  feet,  but  is  usually  about  37  feet.  Two  noise  outputs  are  quoted 
corresponding  to  these  measurements.  The  noise  outputs  for  the  above 
networks  are  shown  in  Table  9. 


TABUS  9 


r-FII/TER  - NOISE  OUTPUT 


Transfer 

Function 

N2  in 

o . 

yd3/seo3 

L - 37  ft. 

L * 50  ft. 

H„(s) 

67.1 

122.5 

^23^S  ) 

46.4 

84.8 

H24(s) 

55.8 

101.9 

Although  the  network  with  transfer  function  H23(s)  gives  better 

results  than  the  other  two , it  does  not  approaoh  the  desired  results . 
With  the  assumed  noise  input,  it  was  necessary  to  allow  the  maximum 
distortion  to  increase  in  order  that  the  noise  output  could  be  brought 
within  bounds.  It  should  be  noted  that  the  maximum  distortion  for 
interoeption  type  courses  occurs  at  short  ranges  where  given  errors  in 
range  rate  cause  smaller  errors  in  lead  angle  than  that  at  long  ranges . 
The  latter  may  be  established  by  examining  the  results  of  reference  (11). 
Although  the  maximum  distortions  for  pursuit  oourses  would  occur  at 
longer  ranges,  these  would  never  attain  the  magnitudes  of  those  for  most 
interoeption  oourses. 
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An  1^3(8)  network,  vhioh  makes  the  mean-square  noise  output  less 
than  25  yda/seoa,  when  L ■ 37  ft.,  is  that  with  transfer  funotion 


Ha3(s) 


s (55  + 45.75s) 

(s  + 2.5)(s  + 4}(s  + 5.5) 


Tabular  results  for  this  network  are  shown  in  Table  10.  The  distortions, 
A0,  were  determined  by  numerioal  integration. 


TABLE  10 

r-FIUTBR  - OUTPUT  ERRORS  - H33(s) 
STRAIGHT  LINES  CASE  (p  = 135°) 


t 

A0  in  yd/seo 

-8 

1.032 

-7.5 

0.943 

-7.0 

0.847 

-6.5 

0.742 

-6.0 

0.625 

-5.5 

0.490 

• -5.0 

0.269 

-4.5 

0.134 

-4.0 

-0.116 

-3.5 

-0.445 

-3.0 

-0.881 

-2.5 

-1.464 

-2.0 

-2.210 

-1.5 

-3.074 

-1.0 

-3.827 

-0.5 

-3.981 

0 

-2.966 

Na  yda/seos 
0 

22.6 

Q 

6.2 

Q 

5.1 
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These  results  are  being  oheoked  by  the  experimental  group,  and  the 
filter  will  be  oheoked  for  distortion  on  other  types  of  signal.  If  it 
is  found  to  be  acceptable  on  other  types  of  signal,  this  network  will 
probably  be  designed  into  the  present  FGCS  breadboard  model. 

It  will  be  necessary  to  oonduot  a more  thorough  investigation  of 
range  rate  filters  when  more  definite  information  concerning  radar 
ranging  noise  characteristics  is  available;  however,  there  would  be 
little  value  in  carrying  the  present  investigation  any  further. 


Bruce  W.  Davis 
Mathematics  Division 


Mathematics  Division 
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NOTATIONS 

A Azimuth  angle  of  the  sight-line  as  measured  from  the  ovnship 

datum  line  i 

8 Slevation  angle  of  the  sight-line  as  measured  from  the  ovnship 

azimuth  plane 

r Present  range  to  the  target 

A,JS,r  Time  derivatives  of  A,  8,  r 
A(t)  The  noise-free  input  signal 

N(t).  The  noise  superimposed  upon  the  input  signal  • 

V_  Ovnship  linear  velocity 

G 

V,j,  Target  linear  velocity 

V Target-ovnship  relative  velocity 

P Angle  betveen  the  target  and  ovnship  linear  spaoe  paths 

A(t)  Time-derivative  of  A(t) 

AQ(t)  Actual  noise-free  output  signal  of  the  filter 
R(t)  Autocorrelation  funotion  of  the  input  noise 

RQ(t)  Autocorrelation  function  of  the  output  noise 

G(f)  Speotral  density  of  the  input  noise 

GQ(f)  Speotral  density  of  the  output  noise 

a Rmpirioal  oonstant  appearing  in  R(t) 

a?  Mean- square  amplitude  of  the  input  noise 

N*  Mean-square  amplitude  of  the  output  noise 

l(t)  The  veighting  or  memory  function  of  a filter 
H(s  ) Laplace  transform  of  w(t) 

A6  Distortion  in  the  output  signal 
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A0y  Maximum  value  of  A6  with  respeot  to  time 

£ 

A0  Mean-square  of  A0  for  0 <.  t < T 

Q Value  of  y/te  * + P 

u o 

Qy  Minimum  value  'of  Q with  respect  to  filter  parameters 

Q Value  of  y/  AO*  +1P 

Z-R  « Zadeh-Ragazzini 

T*  Time  constant  of  the  Z-R  Filter 
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